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. Abstract. Starting from the Maxwell equations, we used the reductive perturbation 

Oh I method to derive a system of two coupled nonlinear Schrodinger (NLS) equations for 

■ the two Beltrami components of the electromagnetic field propagating along a fixed 

, direction in an isotropic nonlinear chiral metamaterial. With single-resonance Lorentz 

1 ^ i ' models for the permittivity and permeability and a Condon model for the chirality 

parameter, in certain spectral regimes, one of the two Beltrami components exhibits a 

I negative-real refractive index when nonlinearity is ignored and the chirality parameter 

QQ ' is sufficiently large. We found that, inside such a spectral regime, there may exist 

T-H , a subregime wherein the system of the NLS equations can be approximated by the 

£J2 [ Manakov system. Bright-bright, dark-dark and dark-bright vector solitons can be 

. ■ formed in that spectral subregime. 
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1. Introduction 

Metamaterials exhibiting a refractive index with a negative real part in certain spectral 
regimes have been a topic of intense research activity over the last decade [T]-|3]. A 
sufficient condition for the exhibition of a negative-real refractive index (NRRI) by an 
isotropic dielectric- magnetic material is that both its permittivity e{uj) and permeability 
jl{uj) have negative real parts for the same value of the angular frequency u [1]. That 
demanding condition can be relaxed if the isotropic material possesses chirality. In this 
case, one of the two refractive indices exhibited by such an isotropic chiral material 
can have a negative real part [3,16]; besides, the possibility of realizing NRRI materials 
through the routes of gyrotropy and nonreciprocity [H Sec. 2.3.4] was considered in 
the more general context of Faraday chiral materials [5]. Although the electromagnetic 
(EM) properties of isotropic chiral materials had been studied in detail in the recent 
past [814T0]. the possibility of isotropic chiral NRRI materials gave fresh impetus [TT] . 
Significant theoretical progress has been reported p!2]4T6] in comparison to experimental 
progress [17]. 

All of the foregoing NRRI materials have been considered to have linear response 
properties. Nonlinear NRRI materials — of the isotropic achiral type — were introduced 
p!8|[T9] shortly after their linear counterparts. Theoretical research has been extensively 
performed on these nonlinear NRRI materials, as they feature nonlinearity- induced 
localization of EM waves and soliton formation [201427] . Similar effects had been 
predicted earlier for nonlinear isotropic chiral materials with positive-real refractive 
indices [28143T] . The possible incorporation of nonlinearity in isotropic chiral NRRI 
materials motivated us to investigate soliton formation in these materials. 

In this paper, we report on the propagation of an EM field along a fixed direction 
in an isotropic chiral NRRI material with nonlinear permittivity and permeability of 
the Kerr type [32]. Furthermore, our analytical approximations may also be useful 
for analyzing field propagation along the distinguished axis of a nonlinear uniaxial 
bianisotropic material. We start from the Maxwell equations in the time domain and 
use the reductive perturbation method (RPM) [33|[M] to derive a system of two coupled 
nonlinear Schrodinger (NLS) equations for the left-handed and right-handed Beltrami 
components of the EM field [9|. Then, we adopt (i) the Lorentz model for the linear 
parts of the relative permittivity and permeability and (ii) the Condon model for the 
(linear) chirality parameter; see, e.g., Ref. [35] and references therein. 

Now, if the chirality parameter is sufficiently large in a certain spectral regime and 
both dissipation and nonlinearity are ignored, the refractive index for the left/right- 
handed Beltrami component is real and negative but that for right /left-handed 
Beltrami component is real and positive [13]. Inside that spectral regime, perturbative 
nonlinearity of the chosen kind induces a spectral subregime wherein the system of the 
NLS equations can be approximated by the Manakov system [36]. Since the latter is 
known to be completely integrable [37H39] , we can predict various classes of exact vector 
soliton solutions that can be supported in the nonlinear isotropic chiral NRRI material: 
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depending on the type of nonlinearity, namely self-focusing or self-defocusing, we find 
that bright-bright solitons (for the focusing case), as well as dark-dark and dark-bright 
solitons (for the defocusing case) can be formed. In all cases, these vector solitons are 
composed of a Beltrami component with negative-real refractive index and another with 
a positive-real refractive index. 

The paper is organized as follows. In Sec. [2] we present the constitutive relations 
and use the RPM to derive, from the time-domain Maxwell equations, the system of 
two coupled NLS equations. In Sec. [3] we show that the coupled NLS equations can be 
approximated by the Manakov system, and we present various types of vector soliton 
solutions. We conclude in Sec. HI 



2. Coupled nonlinear Schrodinger equations 

2.1. Constitutive relations 

In the time domain, the Tellegen constitutive relations of an isotropic chiral material 
are expressed as P-flO] 

D(r, t) = f [e(r) * E(r, t - r) + ^(r) * H(r, t - r)] dr , (1) 

J — oo 

B(r, t) = f [/i(r) * H(r, t - r) - ^(r) * E(r, t - r)] dr , (2) 

J — oo 

where D(r,t) and B(r, t) are the electric flux density and the magnetic induction field, 
respectively; e{t) is the inverse Fourier transform of e(w), and nit) of jliuj)] and ^(t), 
the inverse Fourier transform of ^(o;), accounts for the chirality of the chosen material. 
Hereafter, r = xiL + yy + is the position vector and t denotes time. 

We specify the permittivity and the permeability in the frequency domain in terms 
of linear and nonlinear components as follows: 

~e{uj) =Q(a;) + e„,(|E|MH|2), (3) 
/i(w) = /i^(w) +/i„£(|E|MH|^), (4) 

The linear components of q and jli depend on u and can be decomposed into real and 
and imaginary parts, i.e., 

q(w) = e/j(w) + ze/(w), (5) 
Jj,t,{u) = jxR{u) + ifii{uj), (6) 

where the imaginary parts e/ and jxi{<jj) account for dissipation. The nonlinear parts 
e„£ and jlni are assumed to be independent of frequency and depend quadratically but 
weakly on the magnitudes of both the electric field and the magnetic field. In particular, 
we consider an isotropic Kerr nonlinearity with 

(7) 
(8) 





= ai E 






fi-nt 


= ^2 E 




|H|' 
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where ai^2 and /3i^2 are scalar Kerr coefficients governing the magnitude of the 
nonhnearity. These four coefficients are the necessary ones to describe nonhnear wave 
propagation, in the framework of the NLS model, in an isotropic chiral material; see 
Sec. 2.3 later in this paper along with Ref. [29] for earlier work. The assumed form of 
the nonlinear dependence of Ini and jlni may stem from a Taylor expansion of a more 
general type of nonlinearity (e.g., a saturable one) [13]. 

The frequency-domain chirality parameter ^ (cu) is purely linear, and is defined as 

^{u) = -k{uj) = - [kr^u) + ikj{u)] , (9) 

where c is the speed of light in free space (vacuum), and the frequency-dependent 
function R defines the chiral properties of the material, with its imaginary part kj 
accounting for dissipation together with ij and fij. 
Finally, we define the complex wave numbers 




^±(^) = ^ I /fMAM ± ^(^) I = ~k^^u;) + t~kf{cu), (10) 

and the complex refractive indices 

fi^{uj) = -k^{u) =h^{u) + ihf{uj) , (11) 

CO 

where the superscript + identifies the left-handed Beltrami component and the 
superscript — identified the right-handed Beltrami component, whereas eo and fio are 
the permittivity and permeability of free space. With an exp{—iut) time dependence, if 
E(r,a;) and H(r,a;) are the electric and magnetic field phasors, then the two Beltrami 
components given by 

Q±(r, u) = ^ [E(r, a;)±zr/,(a;)H(r, u)j , (12) 
obey the equations 

V X Q^(r, u) = ±k^{u) Q^(r, w), (13) 
where fjj{u) = fli{u) / ei{u) , and the nonlinear properties have been ignored PITO]. 



2.2. Propagation along the +z axis 

Faraday and Ampere-Maxwell equations are, respectively, stated as V x E(r, t) = 
-{d/dt)B{r,t) and V x H(r,t) = {d/dt)'D{v,t). Let us assume that the EM field 
is propagating along the +z axis. It then has to be transversely polarized: E{z, t) = 
^Ex{z,t) + yEy{z,t) and li{z,t) = iH^^z^t) + yHy{z,t). Furthermore, the field has a 
carrier angular frequency Uc- 

Following the Bohren decomposition of the EM field phasors into left- and right- 
handed Beltrami components [9l|29] provided in Sec. 12.11 we represent 

e-^^^* + c.c., (14) 
e-*'^=* + C.C., (15) 



E{z,t) 



.g+(z,t)e*'==^ + e_ 



q {z,t)e 



H{z, t) = e+p+{z, ty^^ ' + e_p-(z, t)e 
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where "c.c." denotes the complex conjugate, kf = k^iuc), e± = (x ± iy)/v^, 
whereas the field envelopes q^{z,t) and p^{z,t) have to be determined. We also define 

kt'^^^ = dk^{u)/duj\^=uj^ and /c^^^^ = (fk^{u)/duj'^\uj=ujc- 



2.3. Nonlinear evolution equations 



Nonlinear evolution equations for the unknown field envelopes can be found by the 
reductive perturbation method [331IM] as follows. We introduce the slow variables 

Z = e^z, T^ = e{t- kf^^h) , (16) 

where e is a formal small parameter, defined as the relevant temporal spectral width of 
the nonlinear term with respect to the spectral width of the quasi-plane- wave dispersion 
relation [iQl-SS]. Furthermore, we express q'^iz, t) and p'^{z, t) as asymptotic expansions 
in terms of the parameter e as follows: 

g±(Z,T±) = q^{Z,T^)+eqt{Z,T^)+e'qt{Z,T^) + ..., (17) 
p^{Z, T±) = p±(Z, T±) + epfiZ, T±) + e^pt{Z, T±) + . . . . (18) 

Next, we substitute Eqs. (fT7|) and (fT8|) as well as the constitutive relations into the 
Faraday and the Ampere-Maxwell equations, and expand q(ci;), fleiu), and (w) about 
the angular frequency Uc- Assuming that all members of the set {ai^2, (3i,2,^i, 
are of order O(e^), we obtain the following equations at various orders of e: 



0, 



tiW^y drift 



(W^)"9|.fo± 



Here and hereafter, the prime denotes the derivative with respect to uj, 

J e {0,1,2} 



^0 



B 



Tkt + '^fiR 

'-^Ir ±kt - ^KR 



c 



/^o[«2(|g(fP 

eo[ai(|gc)'P ■ 

—iRj jlj 
li iki 



\% 





1 


, J = 




-1 



(19) 
(20) 



(21) 

(22) 

(23) 
(24) 
(25) 



To proceed further, we note that the compatibility conditions required for Eqs. f ll9p - 
( 12T|) to be solvable, known also as Fredholm alternatives [33lll0], are L^W^f^*^]^ 2 = 0; 
where = [l,^iZi] is a left eigenvector of W^, such that L^W*^ = [0,0], with 
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= v^/ifl/efi being the linear impedance, when dissipation is small enough to be 
ignored. Note that this impedance is independent of the chirality parameter h |10] . 

The zeroth-order Eq. ( fT9l) provides the following three results. First, the solution 
f^*" has the form: 

fo± = R±0±(Z,T±), (26) 

where the scalar function (p^i^Z^T^) has to be determined and = [1, =i=zZ^^^]"^ is 
a right eigenvector of W^, i.e., W'^R''^ = [0,0]""". Second, by using the compatibility 
condition L^'^W^'^f^'^ = and Eq. (126!) . we obtain L='=W='=R='= = 0, which is actually 
equivalent to Eq. (ITOl) . Third, the respective zeroth-order Beltrami components of the 
electric and magnetic fields are proportional to each other, namely, = ±iZ£p^ [TO] . 

Next, at 0{e'^), the compatibility conditions for Eq. ([20]) result in L=^(W=^)'R=^ = 0; 
equivalently, 

2{ckf T ujc~^R){ck^ T ^~^r)'\^=^^ = (w'e/?/ii?yL=^,- (27) 

The foregoing equations yield the group speeds = l/kt^^\ which can also be derived 
from Eq. ( |T0]) . Furthermore, Eqs. (120]) and ( 126]) together indicate that ff*" has the form: 

= tiK^)'dT±<p^{Z, T±) + R±^±(Z, T±), (28) 

where the scalar function ip^{Z, T"^) can be determined at a higher-order approximation; 
see, e.g., Ref. [27] . 

Next, at order 0(£^), the compatibility conditions for Eq. ( 12T]) . combined with 
Eqs. (126]) and (128]) . yield the coupled NLS equations 

tdz<P^ - Ikf^'^ dl^<P^ + 7 (101' + \r?) <P^ = -^f ±</)±, (29) 

I /2'\ _____ 

where kc can be obtained from Eq. (127]) . In this equation, the nonlinearity coefficient 
7 = y [eo(aiZ, + AZ^i) + /io(a2^f + /^s^,"')] (30) 

depends linearly on the scalar Kerr coefficients 0:1^2 and /3i^2, and the loss coefficients 
are given as follows: 

f± = c ± ^ V (31) 



2.4- Relation with the Manakov system 

Knowledge of 0^ obtained from solving Eqs. ( 12 9 p immediately yields = 0^ and 
= =pzZ^^0^, by virtue of Eq. (126]) . similarly to the case of a linear material. 
Let us now analyze Eqs. ( 129]) in more detail. First, measuring length z, retarded 
time T = t — z/v'^, and the intensities in units of a characteristic length 

L"*" = tl/lkt^"^^], a characteristic time to, and L+/I7I, respectively, we reduce Eqs. (129]) 
to the following dimensionless form: 

,dz^+ - ^9^0+ + a (|0+|2 + |0-|2) 0+ = -^r+0+, (32) 
z (9z0" - (59t0~) - ^9|0" + a (|0+|2 + |0-|') 0" = -zr-0-. (33) 
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The parameters involved in the foregoing equations are defined as follows: 

u+w L-(2) .+(2) 

I rhC I I r^C I I I 

a = 7/|7|, r± = L+f±. (34) 

Since all nonlinearity coefficients in the ± components of Eq. (129]) are equal to 7, 
the latter parameter was absorbed by our normalization; this way, the nonlinearity 
coefficients in the dimensionless Eqs. fl32|) are fl33|) are equal to the sign of 7. 



Furthermore, as s = sign 



the foregoing NLS system is characterized by (i) 
the parameter 5, which represents the mismatch between the group speeds of the two 
Beltrami components; (ii) the normalized dispersion parameter d for the field 0~; and 
(iii) the linear loss parameters F^. 

In the most general case (5 7^ 0, (i 7^ ±1 and F^ 7^ 0), the system of Eqs. ( 132|) 
and (133|) is not integrable. The same is true even if F+ = F~ = 0; however, for certain 
values of the linear constitutive parameters and the nonlinearity coefficient 7, one can 
find, either analytically or numerically, various types of vector solitons (such as bright- 
bright, dark-dark, or dark-bright ones), as can be gathered from Refs. [^H5]. 

Provided 5 = 0, d = s and F='= = 0, Eqs. f l5^ and flS^ reduce to the Manakov 
system [36]. In our case, the latter takes the generic form of a vector NLS equation, 
which can be expressed as 



s 



iS^u — ^9yU + cr|u|^u = 0, (35) 
where u(Z, T) = [0^(2', T), (p~{Z, T)]^ . The Manakov system is known to be completely 



integrable [37H39] : in fact, it can be integrated by extending the inverse scattering 
transform method that has been used to integrate the scalar NLS equation [161H7] . The 
Manakov system also admits vector A^-soliton solutions [l8] - [53] . 

In the following section, we show that it is possible to find physically relevant 
conditions allowing us to approximate the general system of the NLS Eqs. f l32|) and fl33l) 
to the Manakov system fl35|) and investigate conditions for the existence of various types 
of vector solitons supported by a nonlinear isotropic chiral material. 

3. Vector solitons 

We considered a certain type of isotropic chiral material that conforms to the single- 
resonance Lorentz models for the linear parts of permittivity and permeability, and to 
the Condon model for the chirality parameter (see, e.g., Ref. [55] and references therein) 
as follows: 

ee{u) =eJl- YT^^) ' (36) 

M^) = /^o fl - ... ..^^..... 1 , (37) 
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Figure 1. (Color online) Left panel: Re(n~) (solid red line) and Rc(?i+) (dashed blue 
line) in the spectral regime {Q.2uj^, l.Ga;^). The constitutive parameters used for these 
data are presented at the beginning of Sec. |31 Three distinct spectral regimes exist 
wherein Re(n~) < and Re(ri+) > 0. Right panel: Magnification of left panel in the 
spectral subregime (l.SSwe, 1.47u;e), wherein not only are Re(n^) < and Rc(fi+) > 0, 
but the NLS Eqs. p2|) and (|33p can also be approximated by the Manakov system (j35p . 



Here, Zq = ^y ^o/^o is the intrinsic impedance of free space; {u^,6^}, 
and {cUk, are the resonance angular frequencies and hnewidths of q, fli, and 
K, respectively; while Up and Um characterize the oscillator strengths of respective 
transitions in q and jlf. Finally, a is the rotatory strength of the resonance in the 
Condon model, which measures the degree of chirality. We note that the Lorentz model 
is a widely used fundamental model, describing the dispersive nature of permittivity 
and permeability, also accounting for the presence of a single resonance frequency. On 
the other hand, the Condon model is generally used to represent the dispersive nature of 
chirality; the Condon model in the form of Eq. ( I38ll is such that the real (imaginary) part 
of k{u) is an odd (even) function of frequency [MllSSj, and obeys causality restrictions. 

Next, we used Eqs. ( 136|) . ( 137|) and ( l38l) . as well as the dispersion relations ( ITOl) . 
to determine the refractive indices [cf. Eq. ( ITTj) ] and the parameters involved in 
the NLS equations [cf. Eq. (!34l) ]. All numerical results presented were computed for 
the following values of parameters: Up = O.duje, oJm = O.Swe, = O.ScUe, = O.Gue, 
= 6^ = 6k = 10~^U):, and a = 0.5uj^/c. Note that all resonance angular frequencies 
and linewidths are normalized to the resonance angular frequency of q, which was 
taken to be the largest resonance angular frequency in the constitutive description of 
the chiral material. Other parameter values led to qualitatively similar results. 

For the chosen parameters with a > 0, spectral regimes exist where the refractive 
indices are h~ < and n"*" > if a > Od = 5 x 10~^u^/c. In the top panel of Fig. [1] 
we show an example of the frequency dependences of n~ and n+ for a = 0.5ue/c. In 
fact, if the chirality parameter is sufficiently large, i.e., a > = O.lu^/c, then there 
exists a certain spectral subregime located to the right of the largest resonance angular 
frequency u^, such that n~ < 0, n"*" > 0, and the NLS Eqs. fl32|) and f l33|) can be 
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Figure 2. (Color online) Left panel: Normalized dispersion coefRcients s (dashed 
blue line) and d (solid red line) of the NLS Eqs. dS!])-®, defined in Eqs. in 
the Manakov regime (l.SSwe, 1.47we). Right panel: Normalized loss coefficients, r+ 
(dashed blue line) and (solid red line) in the Manakov regime. Inset shows the 
group-speed-mismatch coefficient 5. 

approximated by the Manakov system ( l35l) . Hence, this spectral regime is referred 
to as the Manakov regime in the remainder of this paper. In the particular case of 
a = 0.5ue/c, the Manakov regime is {1.35ue, lA7ue)] the dependencies of h~ and n"*" on 
u in this regime are illustrated in the bottom panel of Fig. [H Here we should note that 
increase of a above the characteristic value results in the increase of the bandwidth of 
the Manakov regime; for example, when a = 1, the Manakov regime is (1.35cje, 1.73ue). 

If the sign of a is changed, then a role-reversal occurs in that n"*" < and n~ > 
in the Manakov regime. Although we have maintained a > in the remainder of this 
section, the case of a < can be treated analogously. 

In the Manakov regime d ~ s, while 5 ~ and ^ 0. This is confirmed by the 
data presented in Fig. |2l Whereas (i ~ — 1, the maximum deviation between s and d 
(for u = 1.47aJe) is than 10%. Furthermore, F^ are of order of 10"^ while 6 is of order 
of 10~^. Thus, taking into account that the other coefficients in Eqs. (l32l) and (l33l) are 
equal to ibl, we can safely conclude that, in a first approximation, our system can be 
well approximated by the Manakov system ( l35l) . 

Here it should be noted that, following the foregoing arguments, the (small) effects 
of the group-speed mismatch between the two Beltrami components and the linear losses 
can be studied analytically by means of the perturbation theory for solitons; see, e.g., 
Ref . [56] . However, such a study is beyond the scope of this work. In any case, some 
estimations concerning, e.g., the role of losses can already be made, based on the linear, 
dispersionless part of the NLS Eqs. (132|) -( 133|) : one expects that the Beltrami fields 
0^ decay as exp(— F^Z). To further elaborate on this estimation, however, one would 
employ real world parameter values, which are not currently available (since experiments 
on isotropic nonlinear chiral metamaterials considered herein have not been performed 
so far). 

The reduction of Eqs. ( 132|) and ( 133|) to the Manakov system allows us to predict 
different types of exact vector solitons that can propagate in an isotropic chiral NRRI 
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material. In particular, taking into regard that s = d = —1, there are two different 
cases, depending on the sign a of the nonlinearity coefficient 7: 

(a) (T = +1, i.e., the effective Kerr nonlinearity in the Manakov system is self-focusing, 
and 

(b) a = —1, i.e., the effective Kerr nonlinearity is self-defocusing. 
3.1. Self-focusing nonlinearity in the Manakov regime 

When s = (i=— l,cr = +l, and the boundary conditions associated with Eq. (l35l) are 
0^ — )■ as T — )■ 0, there exist exact bright-bright soliton solutions, i.e., both Beltrami 
components take the form of bright solitons. Following Refs. [49l[50], we may find such 
solutions by using the traveling-wave ansatz 

(/.+ (Z,T) = ni(Oexp [i{l + C"^ /2)Z + iCT] , (39) 
0-(Z,T) = M2(0exp [i{q^ + C''/2)Z + iCT] , (40) 

where ^ = a/2 (T — CZ), the parameter C sets the (common) speed of both components 
of the traveling wave in the (Z, T)-plane, q is an arbitrary parameter connected with 
the wave number of the right-handed Beltrami component, and Ui^2 are unknown real 
functions. 

Introducing the ansatz comprising Eqs. fl39|) and fHOj) into the Manakov system 
( l35|) . we obtain the following system of ordinary differential equations (ODEs) for Ui 
and U2: 

d^ui 



di 

d^U2 



2 -ui + {u{ + ul)ui = 0, (41) 
Y - q^U2 + {ul + ul)u2 = 0. (42) 



d^ 

The simplest exact solution to this system has either ui = or U2 = 0, i.e., either 

Ul = 0, U2 = V2qsech.{q^) (43) 

or 

Ul = V2sech(0, U2 = 0. (44) 

This solution can be generalized to the case wherein both ui and U2 are nontrivial 
[49ll50] . Particularly, if g = 1, a simple one-parameter family of such solutions is 
composed of symmetric, single-humped bright solitons of the following form: 

Ul = -\/2 cos 6'sech(^), U2 = -\/2sin^sech(^), (45) 

where 6 is an arbitrary parameter. On the other hand, if < g < 1, the soliton is 

Ul = v'2(W)cosh(gO 

cosh(A,^) cosh(g^) — gsinh(A^) sinh(g^)' 

= -9^2(1 -g^)sinh(AO ^ ^^^^ 

cosh(A^) cosh(g^) — gsinh(A^) sinh(g^)' 
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where = ^ — £^nd is an arbitrary constant. The solutions fl46p and ( 14 7p are 
generally asymmetric; however, ui is symmetric but U2 is antisymmetric when = 0. 
Other solution branches, such as soliton bound states and wave- and daughter-wave 
solutions (for Ui -C U2 or U2 <C Mi), have also been found; see, e.g., Refs. [l9l[50], as well 
as the book |13] and references therein. 



3.2. Self-defocusing nonlinearity in the Manakov regime 

Let us now change the sign of a to a = —1 and keep s = d = —1. Now, there exist both 
dark-dark solitons (where both the Beltrami components take the form of dark solitons) 
and dark-bright ones (where one Beltrami component is a dark soliton and the other is 
a bright soliton). 

Dark-dark soliton solutions of the Manakov system fl35|) . with boundary conditions 
|(/)^| — )■ /i-t as T — 7- 0, where ^± are the amplitudes of the continuous- wave (cw) 
backgrounds carrying the dark solitons, can be expressed in the following general 
form [5T1I521: 



Z, T) = fi±{cos ip± tanh[i/(T — VZ)] + i sin (p±} 



xexp z ( ^ + + ) Z + zr/±T , (4^ 

where 1/ and V are the (common) inverse width and speed of dark solitons, while rj^ 
are constants connected to the background angular frequency Uc and wave numbers kf. 
The phase angles (p± = tan~^[(?7-|- — V)/^] are connected with the rest of the soliton 
parameters according to the relations z/^ = /i^ cos^ + cos^ while the soliton 
intensities are as follows: 

= l^l-uhechy{T -VZ)]. (49) 

Finally, for the same case (cr = —1), we present the exact analytical dark-bright 
soliton solutions of the Manakov system ( l35l) [^[52]: 

(j)~^{Z,T) = 0o{cosv9tanh[D(T — CgZ)] + icosip} 

X exp[inT + i{n'^ + (l)l)Z], (50) 
0"(Z,T) = r]sech[D{T -CsZ)] 

X exp{2C,r + «[(D2_c2)/2-02]^}, (51) 

where 0"'' and correspond, respectively, to the dark and bright soliton components; (po 
is the amplitude of the cw background carrying the dark soliton; if = ta.n~'^[{Q — Cs)/D] 
is the dark soliton's phase angle; D and Cg denote the (common) inverse width and 
speed of both components; Q sets the frequency of the dark-soliton component; and 
rj is the bright-soliton amplitude, with 77^ = (p^cos'^ip — . In the limiting case of 
?7 — 0, the bright soliton component vanishes and the inverse width of the dark soliton 
becomes D = cpQCosip; this way, Eq. (!50!) describes the single dark soliton solution of 
the defocusing NLS equation (HJHS] . It is important to stress here that the dark-bright 
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soliton of Eqs. (ISO]) and (IHT]) is a generic example of a symbiotic soliton: the bright- 
soliton component exists only due to the coupling with the dark soliton and would 
not be supported otherwise, due to the fact that the NLS equation for 0" features a 
self-defocusing nonlinearity. 

Finally, we note in passing that other soliton solutions of the Manakov Eqs. (l35ll 
can be systematically constructed by employing the symmetries of the system (such as 
the SU(2) rotational symmetry) [53] . 

4. Conclusions 

In conclusion, we have studied the propagation of electromagnetic pulses in isotropic 
chiral NRRI materials with Kerr nonlinearity. We have employed the reductive 
perturbation theory to derive a system of two coupled nonlinear Schrodinger (NLS) 
equations for the envelopes of the two Beltrami components of the EM field. Assuming a 
single-resonance double Lorentz model for the effective permittivity and permeability of 
the medium, as well as a Condon model for the chirality parameter, we have shown that 
if the characteristic chirality parameter is sufficiently large, then the material exhibits 
a negative- real refractive index for the right /left-handed Beltrami component in certain 
spectral regimes whereas the left /right-handed Beltrami component does not. 

Furthermore, we have shown that in a certain subregime inside such a spectral 
regime, the system of the NLS equations (which is non-integrable in general) may be 
approximated well by the completely integrable Manakov system. This way, we have 
predicted various types of vector solitons, such as bright-bright, dark-dark or dark-bright 
ones, that can be formed in the chosen type of chiral material. Let us emphasize that, 
in all cases, the presented vector solitons are composed by a Beltrami component that 
features negative-real refractive index and another Beltrami component that features a 
positive-real refractive index. 

The propagation of the predicted vector soliton solutions needs to be examined by 
direct numerical simulations in the framework of the Maxwell equations. Furthermore, 
a relevant future challenge is the analysis of pulse propagation in the bianisotropic [5ll57] 
counterparts of the isotropic materials we have discussed here. For that purpose, our 
work can be emulated in three steps as follows. In the first step, one has to obtain 
explicit plane wave solutions and associated dispersion relations, similar to the ones in 
Eqs. ( ITO|) . In the secon step, one has to consider propagation of wavepackets, similar 
to Eqs. (IHj) and (fTSl) . along a fixed direction (identified by the z-axis in our work) of 
the linearized material — see, e.g., Ref. [57]. The third and final step is to apply the 
methodology of Section 2.3 to derive a vector NLS-type model (and, in the best-case 
scenario, a Manakov system) describing nonlinear pulse propagation in the nonlinear 
material. Also, another interesting possibility would be to apply the derived results in 
the case of pulse propagation inside a layered chiral [58l|59] or a layered bianisotropic 
material [60l|6T]. For this purpose, one should first find expressions for the EM field 
components inside each homogeneous layer (relying, e.g., on plane waves for the linear 



Vector solitons in nonlinear isotropic chiral metaniaterials 



13 



layers, and nonlinear periodic waves or solitons for the nonlinear layers) and then impose 
the boundary conditions at the interfaces. Finally, our work may be exploited for the 
analysis of pulse propagation inside composite materials comprising linear chiral and 
nonlinear dielectric materials [62]. Such studies are currently in process and relevant 
results will be reported in the future. 
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